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Abstract 

An analysis is made of various methods of phenomenological renor- 
malization based on finite-size scaling equations for inverse correlation 
lengths, the singular part of the free energy density, and their deriva- 
tives. The analysis is made using two-dimensional Ising and Potts lat- 
tices and the three-dimensional Ising model. Variants of equations for 
the phenomenological renormalization group are obtained which ensure 
more rapid convergence than the conventionally used Nightingale phe- 
nomenological renormalization scheme. An estimate is obtained for the 
critical finite-size scaling amplitude of the internal energy in the three- 
dimensional Ising model. It is shown that the two-dimensional Ising and 
Potts models contain no finite-size corrections to the internal energy so 
that the positions of the critical points for these models can be determined 
exactly from solutions for strips of finite width. It is also found that for 
the two-dimensional Ising model the scaling finite-size equation for the 
derivative of the inverse correlation length with respect to temperature 
gives the exact value of the thermal critical exponent. 

PACS: 05.10.Cc, 64.60.Ak, 05.50.+q, 75.10.Hk 
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1 Introduction 



As we know [jy, the renormalization group method lies at the basis of modern 
phase transition theory. In this approach, the critical points of a system are 
identified with the fixed points of renormalization group mappings. By lineariz- 
ing the renormalization group transformations near a fixed point it is possible 
to find the critical exponents of the system of interest. 

Many different specific realizations of the renormalization group procedure 
now exist. One of these realizations is the so-called phenomenological renormal- 
ization group proposed by Nightingale || || (see also the reviews 0, D). In this 
approach the scaling functional relationship between the correlation lengths of 
partially finite subsystems is interpreted as a renormalization group transfor- 
mation. Combining this with expressions for the correlation lengths in terms of 
eigenvalues of the transfer matrices reduces the initial functional equation to a 
transcendental one and its solution gives an estimate for the critical tempera- 
ture. 

Dos Santos and Sneddon Q noted that a phenomenological renormalization 
group can be constructed using not only the correlation length but any quantity 
having a power divergence at the phase transition point. 

Binder [Q, |) suggested a variant of the phenomenological renormalization 
group based on using the moments of the distribution function of the order 
parameter; the moments can be expressed directly in terms of the susceptibility 
of the system. This method has been widely used in calculations using the Monte 
Carlo method. Recently Itakura Q expanded the Binder phenomenological 
renormalization scheme using linear combinations of different-order moments. 

In the present study various phenomenological renormalization variants us- 
ing a wide range of physical quantities (such as the free and internal energy, 
ordinary and nonlinear susceptibility, and so on) are used to calculate the crit- 
ical coupling in the two-dimensional Ising and Potts models and also in the 
three-dimensional Ising model. An analysis can reveal the best renormalization 
strategies. Increasing the convergence is fundamentally important for three- 
dimensional systems because as a result of the giant orders of the transfer matri- 
ces, the eigenvalue problem can only be solved for subsystems having extremely 
small transverse dimensions. 

However, a study of the convergence of various types of renormalizations 
revealed equations which can be used to find exact values of the parameters of 
infinite systems (critical temperatures and in one case, critical exponent) using 
the characteristics of partially finite subsystems. 

This article is constructed as follows. Section 2 contains equations for differ- 
ent variants of the phenomenological renormalization group. In Section 3 these 
equations are applied to the two-dimensional Ising and Potts models and in Sec- 
tion 4 they are applied to the three-dimensional Ising lattice. The conclusions 
reached are summarized briefly in the final section. 
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2 Equations for the phenomenological renormal- 
ization group 

In accordance with finite-size scaling theory |R| [□], g [HJ, the inverse corre- 
lation length kl and the singular part of the reduced free energy density /£ 
{fl = II — foo, , where /l and /oo are the free energy densities of the sub- 
system and the total system, respectively) near the critical point t = h = 
(t = K — K c is the "temperature" or more accurately the deviation from the 
critical coupling, and h is the external field) satisfy the functional relationships 

K L {t,h) =b- l K L/b (t\ti), (1) 

fi(t,h) = b- d fl /b (t>,h>). (2) 

Here d is the space dimension, L is the characteristic linear size of the subsystem, 
b = L/L' is the coefficient of the scaling transformation. The renormalization 
group equations t' = t'(t,h) and bl = h'(t,h) linearized near the fixed point 
t = h = have the form 

t' = b yt t + 0(t 2 ), (3) 

ti = b Vh h + 0{h 2 ) 1 (4) 

where yt and yu are the thermal and magnetic critical exponents of the system. 

In the Nightingale approach |^|, ||, |J, |) the relationship ([!]) is treated as 
a renormalization group mapping (t,h) — > (t',h'). In this case, the inverse 
correlation length of the spin-spin correlation function is given by the standard 
expression 

K L = \n(\[ L) /\\i L) \), (5) 

where \\ L ^ and are the largest and second largest eigenvalues of the transfer 
matrix of the subsystem, respectively. As a result, for systems with phase 
diagram symmetry h — > —h (and thus h c = 0; we only consider these systems 
here) the only unknown coordinate of the critical point can be determined from 
the transcendental equation 

Lk l (K c ) = (L-1)k l ^(K c ), (6) 

(in order to maximize the accuracy of the estimate K c it is usually assumed that 
L' = L — l). Nightingale described his renormalization method as "phenomeno- 
logical" H (see also 1 14 ) since it uses a scaling finite-size equation which is not 



derived microscopically using the method. 

Quite clearly, the phenomenological renormalization (t, h) — ► (t! ,hf) can be 
made with equal justification using relation (^). In this case, we obtain the 
following equation to determine the critical coupling K c 

L d fl(K c ) = (L-l) d fl_ 1 (K c ). (7) 
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The free energy density of the partially finite subsystem L d 1 x oo is then 
calculated using the formula 

f L = L 1 - d \n\[ L \ (8) 

and the "background" /oo is taken from outside. 

Differentiating expressions (1) and (2) m times with respect to temperature 
and n times with respect to the held allowing for (3) and (4), we obtain [15] 

L l-my t -ny hK (rn,n)^ ^ = K (m,n) ^ ^ ^ 

L d-myt-ny h ^s[m,n) ^ ^ = ^ d -my t -ny h jsjm,n) ^ yy (10) 

These relationships or some combinations thereof can also be considered as 
implicit renormalization group transformations. For large subsystem sizes all 
equations of this type should give the same results. However, for small L which 
we need to deal with in practice, the various renormalization group equations 
yield estimates of different accuracy. As we shall see, this circumstance means 
that it is possible to construct phenomenological renormalization schemes which 
give faster convergence than that achieved in conventionally used realizations 
of the phenomenological renormalization group method. 

In addition to equations (||) and (Q), in the present study we also give esti- 
mates of the critical coupling obtained from the following renormalization group 
equations [taking those ratios and jl0| ) , from which the unknown exponents 
yt and yh are dropped a priori] : 



(4) 



L d x 



(4) 



-1 



K c (L - lyxl 



(11) 



where xl = 9 2 /l/9/i 2 |/ 1= o = fj}° is the susceptibility and Xx = d 4 f l I dh A \h = $ 

is the nonlinear susceptibility (equation ( pi] ) corresponds to the Binder 
phenomenological renormalization |7|, |8J); 

L"- d {^flxL = {L— if-o^f/xL-i; (12) 

l'-'^/xl = (L — ly-^ljxL-n (13) 

L^nf/xl = {L-lf-^f_J X l_ 1 . (14) 

Here we have — d n KL/dh n \h=o — ft^'™' (K, 0). Formulas for the derivatives 
of the inverse correlation length and the free energy with respect to h, expressed 
in terms of the eigenvalues and eigenvectors of the transfer matrices, suitable 
for programming are given in [ fl6| [l7| . 

For better clarity equations (gh (@), @), (ph, @, and (p]) will be de- 
noted by the mnemonic symbols "k", "/ s ", "jc^/x i u {^ 1 ') /X, " k( * 2 Vx'\ 
and "k*- 4 ) /x 2 " , respectively when the numerical data are given in the tables. 
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3 Two-dimensional systems 



We shall first test various types of phcnomcnological renormalization group 
equations on the two-dimensional Ising and Potts models, for which the positions 
of the critical points are known in exact analytical form 18, [l9|. 



3.1 Ising model 

We shall start with the most studied system, the Ising model. Table [l] gives 
results of our calculations of the critical coupling for a simple square Ising lat- 
tice. The calculations were made for L x oo strips having a periodic boundary 
condition in the transverse direction. Estimates are given for (L — 1,L) pairs 
with L < 5. For a (3,4) pair the relative errors of the estimates are also given 
in parentheses; their sign indicates whether they are lower (-) or upper (+) esti- 
mates. The types of renormalization group equations used are indicated in the 
first column. We recall that in this model we have 

i oc 1 

K c = -ln(l + 72), f x = — + -ln2 

Z 7T A 

(G is the Catalan constant) [20] . 

It can be seen from the data presented in Table [j] that the best lower estimate 
of K c is obtained from equation (|lj) . The Binder scheme gives slightly inferior 
results. This approach, which is usually used for Monte Carlo simulations, as we 
noted in the Introduction [because it does not have formulas of the type (||) for 
the inverse correlation length], was used in the transfer matrix variant in plf . 
The Nightingale renormalization (first row in Table [j]) which is almost universal 
(except perhaps for the study just mentioned pi)]) is used for calculations using 
the transfer matrix method and is only third in terms of accuracy among the 
lower estimates. 

It is worth noting that we also succeeded in obtaining phcnomcnological 
renormalization group equations which also yield upper estimates of K c (the 
two lower rows in Table [l]). Of these two estimates the renormalization scheme 
using the singular part of the free energy density [equation(^)] gives better- 
quality results. The values obtained using this scheme are at the same time the 
best in terms of absolute value among all the lower and upper estimates. 

Unfortunately the renormalization procedure using the singular part of the 
free energy requires an advance knowledge of the critical free energy of the 
system /oo- In order to overcome this constraint, instead of the two strips 
required to construct the renormalization group equation, we take three. In 
other words we shall take the triads (Li, L2, £3); an additional relationship 
can eliminate the parameter /oo- For the free energy densities of the three 
subsystems at the critical point we have 

f Ll =foo+A f /Li, f La =f 00 +A f /Li f La =f 00 +A f /Li (15) 
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Here Af = L /£(0, 0) is the critical finite-size scaling amplitude of the singular 
part of the free energy density. Eliminating the parameters foo and Af from 
the system (fLo), we arrive at the equation 



L d 2 -Li(L 



L3 — ^2 



where 



F(K) = { Ll f f\ (17) 

JL 2 ~ J L 3 

Equation ( |l6| ) can apparently be used to determine the position of the critical 
point but it has no real roots. 

All the eigenvalues of the transfer matrix of an L x 00 Ising cylinder are 
known in exact analytic form for any L [po[ p2[ . The largest eigenvalue is 

\[ L) = (2sinh2^) L / 2 exp[( 7l + 7 3 + .--+72L-i)/2], (18) 

where j r are the positive solutions of the equations 

cosh 7,. = cosh2_ftTcoth2i<" — cos(irr/L). (19) 

Using formula (|8|) we then obtain expressions for the free energy densities 
(i = 1, 2, 3) contained in (p7|). 

Subsequent calculations using these formulas show that the curves T{K) do 
not intersect the abscissa and thus the equation T(K) — has no real solutions 
(see Fig. 1). However, it is noticeable that for various choices of parameters 
Li, L2, and L3 the function F(K) [and therefore F{K)\ has an extremum 
positioned at the exact value of the critical coupling. Differentiating relation 
( [l7|) with respect to if, we arrive at the equation 

(u Ll - UL 2 )(/L 2 - h 3 ) = (UL 2 ~ U L3 )(f Ll - /l 2 ), (20) 

where ul = dfL/dK is the internal energy of the subsystem. By using this 
equation or finding the extremum of the function F(K) we can extract exact 
values of the critical temperatures for infinite systems (at any rate for some) 
from the characteristics of their finite subsystems. (This is the desired aim of 
any finite-size theory.) 

Equation (20|) can be approached from several positions. We shall analyze 
expression (|l0) in the first derivative with respect to K (m = 1, n = 0). For a 
fixed point this rcnormalization group transformation gives 

L d -y<u s L (K c ) = (L') d -y*ui,(K c ), (21) 

where u s L = ul — Uqo is the singular part of the internal energy density. It 
follows from ( |l9| ) that 

d lr /dK\ K=Kc = (r^O) (22) 
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(in this context see Fig. 118 in the book ^3|). After differentiating we find that 
ul(K c ) = \/2 regardless of L. In other words, all the finite-size corrections to 
the background Moo (= V2 p(J) are zero, i.e., the singular part of the internal 
energy is completely absent. Consequently, we obtain 

u' L (K B ) = 0. (23) 

This property of the model leads to equations ( |20| ) and (|2l]) being satisfied. 

We shall now analyze the first derivative with respect to temperature of the 
other fundamental equation in finite-size scaling theory, i.e., equation (1). 

The second eigenvalue of the transfer matrix of the Ising model in terms of 
value is given by 

A^ L) = (2sinh2^) L / 2 exp[( 70 + 72 + ... + 72L-2)/2]. (24) 

Equations (5), (18), (19), and (24) can be used to find an analytic expression 
for the inverse correlation length. Taking the first derivative with respect to 
K we obtain the dependencies Kl(K) = Okl/OK for various L. All these 
dependencies have a common point of self-intersection and this is positioned at 
the exact value of the critical coupling: k^iKc) = —2 for all L — 1,2,... (i.e. 
in an infinite lattice koo = —2). Thus we have 

k L {K c ) = k L ,(K c ). (25) 

Reversing this reasoning, we conclude that equation (|25|) can also be used to 
determine the exact position of the critical point of an infinite two-dimensional 
Ising lattice using the solutions for strips of finite width. 

Moreover, in accordance with relation (^) for m — 1 and n — we have 

Consequently, this finite-size expression taking into account condition ( p5| ) gives 
an exact value for the thermal critical exponent of the two-dimensional Ising 
model: y t = 1. 

For the Ising model equation (|25| ) can be generalized to the form 

(k Li - k L2 )(K L2 - k Ls ) = (k L2 - k L;s )(n Ll - k L2 ). (27) 

and the problem of determining K c can be reduced to searching for the ex- 
tremum of the function 

$ = KLl 1 KL \ (28) 

K L 2 K L 3 
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3.2 Potts model 



In the two-dimensional q-state Potts model for an isotropic square-cell lattice 
the critical coupling is K c — ln(l + ^/q) pi| . Table || contains estimates of K c 
obtained by us using various constructions of the phenomcnological renormal- 
ization group for the Potts model with three states (q — 3). The free energy 
background in this model is given by 

/«, =4G/37r + ln(2V3) + iln(2 + V3). (29) 

An analysis of the data given in Table || shows that equations (^) and ( ft2| ) 
give substantially better-quality estimates than the Nightingale approach (first 
row in the table). The renormalization group equation based on the scaling 
relationship for the singular part of the free energy again gives the smallest 
errors in terms of absolute value. 

Calculations made for Lxoo strips with a periodic boundary condition in the 
transverse direction then uniquely indicate that in the two-dimensional q-state 
Potts model all finite-size corrections to the internal energy (uqo = 1 + 1/yfq 
[jisl ) are also absent at the critical point. Consequently, the equation 

u L {K c )=u L ,{K c ) (30) 

gives the exact value of K c . Satisfying condition ( J30| ) then has the result that 
equation ( |2C| ) is valid and thus the extremum of the function ( |l7| ) is positioned 
strictly at the exact value of the critical coupling. It should be noted that in 
computer calculations it is more logical to search directly for the extremum of 
the function than to differentiate the free energy numerically and then solve the 
transcendental equation ( |30| ) numerically. 

Equation ( |30j ) for the pair (1, 2) was first obtained by Wosiek from his postu- 
lated maximum principle for the normalized moments of transfer matrices p7| 
(see also p8[). In the Wosiek equation was then generalized to the form 
(|30|). Our equation (20) is a further generalization of (|30|). 

Note that in the Potts model the curves kl(K) also intersect. However, 
unlike the Ising model the points of intersection are not now positioned at the 
exact value of K c . At the exact value of the critical coupling the functions 
kjiK) have different values for various L. Unfortunately the scaling formula 
(|26|) then only gives approximate estimates of the thermal critical exponent 
which, for example for the three-site Potts model is yt = 6/5. 

Returning to the two-dimensional Ising model we can confirm that a com- 
plete system of equations of the type ( [sof ) is in fact obtained for this model: 

ui(K c ) = ui,(K c ), (31) 

where u 3 L = df J L /dK are internal energy levels; this term was taken by analogy 
with the "free energy levels" f[ = L 1 ^ d ln\^ L \ j = 1,2,... (see, for example 
[30]). 
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4 Three-dimensional Ising lattice 



We shall now consider a three-dimensional Ising model for a simple cubic lattice. 
Unlike the two-dimensional Ising model, which is solved analytically over the 
entire temperature range, or the two-dimensional Potts model which is solved 
exactly at a single point, the critical point, for the three-dimensional Ising model 
we are forced to confine ourselves to numerical estimates. The highest-quality 
estimates for the critical coupling and the free energy background have now been 
obtained by Monte Carlo simulation for L x L x L cubes: K c = 0.22165459(10) 



Table |3j gives our results for the three-dimensional Ising model. The renor- 
malizations were made for LxLxoo parallelepipeds with periodic boundary 
conditions in both transverse directions. The calculations were made for sub- 
systems with side lengths L < 4, i.e. for transfer matrices having dimensions 
up to 65 536 x 65 536. 

We can see from Table | that as in the two-dimensional case, the best lower 
estimates of K c are provided by the renormalization group equation (|l4|). On 
comparing the lower and upper estimates we can conclude that the best of all 
the results is again given by renormalizing the singular part of the free energy. 

In the three-dimensional case the amplitudes of the finite-size corrections to 
the internal energy at the critical point are no longer zero. In accordance with 
the scaling equation (|2l|) we have 



where A u is the amplitude of the leading finite-size term in the internal energy 
density of the system at the critical point. Calculations using formula (^) for a 
4x4xoo cluster using the values Uoo = 0.990 637(47) fHH andy* = 1.5865(14) 
H| give A u ~ 1.19. 

Since A u ^ 0, equation (|3p| ) now only reproduces an approximate value of 
K c . Also equation (|2(]) is not satisfied exactly. Calculations of the position of 
the minimum of the function ( |l7j) for the triad (2, 3, 4) yield the estimate K c = 
0.233003. Unfortunately, the accuracy of this estimate (5.1%) is significantly 
inferior to the errors of the estimates obtained from the other renormalization 
schemes (see Table ||) . 

5 Conclusions 

Nontrivial models with known values of the critical coupling have been used 
to study phcnomcnological renormalization schemes which for the same sub- 
system sizes give more accurate results than those used so far. We established 
that for Ising models the renormalization group equations using combinations 
£1-2(2^.(4)^,2 anc [ 2y 1_d K^ 2 - ) /xL, give uniformly converging bilateral estimates 



[0 and /, 




0.77790(2) HJ. 



u L (K c ) = Uoc + A u L- d +y* + . . . . 



(32) 
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of K c where the errors (in the three-dimensional, most important case) are less 
than half those obtained in the Nightingale method. 

Another result of this study is related to the exact results. We have shown 
that by making the amplitudes of the finite-size corrections to the internal en- 
ergy vanish, exact values of the critical temperatures can be determined for 
the infinite two-dimensional Ising and Potts models from the solutions for their 
partially finite subsystems. In addition, we have observed that for the two- 
dimensional Ising model the temperature derivative of the inverse correlation 
length has absolutely no finite-size corrections. This means that it is also possi- 
ble to find the exact value of the thermal critical exponent for this model from 
the finite-size equation. 

Both the conclusions reached from these results have practical applications 
but separate analysis is required to develop them. 
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Table 1: Estimates of K c for a two-dimensional square Ising lattice; K% 



0.440 686... 


eq. 


(2,3) 


(3,4) 


(4,5) 


(«) 


0.42236 


0.43088 (-2.23%) 


0.43595 


(x (4) /x 2 ) 


0.42593 


0.43242 (-1.88%) 


0.43672 


(« (4) /x 2 ) 


0.42596 


0.43243 (-1.87%) 


0.43673 


(/ s ) 


0.44324 


0.44168 (+0.23%) 


0.44105 


(« (2) /x) 


0.47420 


0.45153 (+2.64%) 


0.44626 



Table 2: Estimates of K c for a two-dimensional square three-site Potts lattice; 



K exact 1.005Q52. 



cq. 


(2,3) 


(3,4) 


(4,5) 


(«) 


0.96248 


0.98350 (-2.1%) 


0.99467 




0.99311 


0.99920 (-0.6%) 


1.00380 


(D 


1.00927 


1.00667 (+0.2%) 


1.00565 



Table 3: Estimates of K c for a three-dimensional simple cubic Ising lattice; 



K exact 0.221654 59(10) 



eq. 


(2,3) 


(3,4) 


(«) 


0.21340 


0.21826 (-1.53%) 


(x (4) /x 2 ) 


0.21823 


0.22002 (-0.74%) 


(kW/x 2 ) 


0.21824 


0.22006 (-0.72%) 




0.22354 


0.22236 (+0.32%) 


(« (2) /x) 


0.22658 


0.22323 (+0.71%) 
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Figure 1: Behavior of the function T(K) for triads of Ising strip clusters 
(1,2,3), (2,3,4), and (3,4,5) given by curves 1, 2, and 3, respectively. In 
all cases the minimum of the function is positioned at the exact value K c = 
\ ln(l + s/2) = 0.440686 . . .. 
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